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Abstract 

The star-matrix models are difficuh to solve due to the multiple powers of the 
Vandermonde determinants in the partition function. We apply to these models a 
modified Q-matrix approach and we get results consistent with those obtained by 
other methods. As examples we study the inhomogenous gaussian model on Bethe 
tree and matrix q-Potts-like model. For the last model in the special cases q = 2 and 
g = 3, we write down explicit formulas which determinate the critical behaviour of 
the system. For q = 2 we argue that the critical behaviour is indeed that of the Ising 
model on the ch^ lattice. 
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1 Introduction 



The star-matrix models were considered first in j^] in connection with the q-Potts 
model and percolation problem. Another field of interest is the "induced" QCD which 
in the large limit is equivalent with the star-model on Bethe tree. The star-models 
are also a direct generalization of the c = 1 matrix model, which is the particular 
case for g = 1. 

The difficulties in solving exactly these models are related with the multiple pow- 
ers of the Vandermonde determinants in the partition function. In the Q-matrix ap- 
proach, another problem related with the first , is the choosing of a proper basis in 
which to define consistently the Q-matrices. 

In this work we show that it is possible to apply a modified Q-matrix approach 
and to define consistently the Q-matrices. The method gives consistent results for 
the cases where other methods can be applied :the saddle-point method, Schwinger- 
Dyson approach, gaussian integration etc. The gaussian models and matrix model 
on Bethe tree give similar results with the previous ones. However, our method can 
be applied in all genera, and permits a more precise study of particular parts of the 
Bethe tree. 

In the last section we consider the matrix formulation of the q-Potts-like model 
and study it as a polymer on a random surface. It is important that we still get 
consistent results for this model. It is a non-trivial result, because the coupling 
conditions are overdetermined: the Q-matrices depend on 3g + 3 variables and we 
have 4(7 + 3 equations. More generally when the central potential Vq is gaussian and 
the lateral potentials Va are of order n we have 3q + n variables and (n + l)g -|- 3 
equations always with combatible solutions. 



2 Q-matrix approach 

The Q-matrix approach was introduced for the 1-matrix model and further de- 
veloped for the two-matrix and the chain multi-matrix models in Here we are 
applying this approach to another class of multi-matrix models - the star-matrix 
models. 

The partition function of the g-star model is given by: 

„ <? q q 

Z= Y[ dM^dMo exp(Vb + Y.^a + MoJ2 '^M (2-1) 

■' a=l OL a=\ 

with the potentials Va = J2r=i ^a,rM^,a = 0,1 ... q. It is possible to integrate over 
the angular degrees of freedom and to remain with integrals only over the eigenvalues: 

I- N q q 

Z= X{{d\o,i n <iA,,i)(A(Ao))'"'^ n ^i^o^V (2.2) 

with 

N q q 

y = E(^o(Ao,i) + 

1=1 a=l a=l 
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We define the orthogonal functions basis as and q + 1-conjugate functions ria,m'- 

Sm,m^ (2.3) 

a=l a=l a=l 

The basic functions have the property: 

^°'"7?o,m(Ao) = / dA„e^"+^"^o^"77„,^(A„) (2.4) 



Introducing this relation ( |2.4| ) in the orthogonality condition (2.3) we get: 



J dAoeo'n(Ao)e^« n Vlmj^o) =hnf[ Sn,n.^ (2.5) 

Because the basic functions Co,n(Ao), f/o,n(Ao) at power q are linear combinations of 
basic functions: 

nq nq 

^0,n(Ao) = ?0,n(5r(Ao) + ^ afc^O,nij(Ao), ^o,n(Ao) = VO,nq{Xo) + ^ CLk'r]0,nq (Ao) (2.6) 
k=l k=l 



we can show that the integral ( |2.5D follows from the orthogonal condition : 

J C?Ao^o,n(Ao)e^''r/o,m(Ao) = K5n,m (2.7) 



and: 

nq 



hn — no^nq + ^ | Ofc | ^/lO,nq-fc 



k=l 



Inserting the expression of r/o,m ( |2.4D in the relation (^^)we get the orthogonal con- 
dition: 

" dXodX^^o,n{Xo)e^°^^''+'-^''^-v'a,miK) = K,nSn,m (2.8) 



The property (^^J) gives the possibility to integrate over Vandermonde determi- 
nants: 

N-l , . 

n T^det[%,.(Ao,,)] = / dA„e^-+'^"^»^"det[r?„,,(A„j)] 
and permits to calculate the partition function: 

N-l I q \ 

Z = constm n hl-^ n K^, (2.9) 

i=0 ' \a=l / 

We introduce the Q-matrices as: 

/ dXo n dA„C^(Ao)A„e^ [] 

(2.10) 

a=l o=l a=l 



and the P-matrices as: 

d 



f dXo n dXa^'H^o)^ (Cn(Ao)e^^+"«S'^i) eSLi^'^ f[ r/„,^JA«) = (2.11) 

a=l a=l 

q 



Q = l 



a=l a=l 1^ a^Q,f5 



q 

J 



We can now derive the couphng conditions: 



gPo + K(Qo) + E CaQa = (2.12) 
a=l 

Pa + VaiQa) +CaQo = 0, a = I, . . . q 

We consider only the symmetric case when the order of potentials is pa = P- The 
calculation of degree for matrices gives: 

Qa G [— "T-i^^li for all a 

Qo G [-mo, no], (2.13) 

where m = l,n = po — 1; mo = p — 1, no = 1. 

The notation shows that we have Qo-matrices with finite band with m lower and 
n higher diagonals. 

We consider only 1-point correlation functions, hence we do not need the flow 
equations. The 1-point correlation functions can be calculated in the usual way as: 

< TtM^^ . . . M^" >= T¥(Q^i ■■■Qg") (2.14) 

These 1-point correlation functions can be calculated in every genus h and in principle 
for arbitrary potentials, not only gaussian. 



3 The gaussian model 

After the exact solution of c = 1 matrix model (or chain q-multimatrix model) 
with gaussian potential, the star q-multimatrix model is an obvious target for study. 
Even with a trivial gaussian potential, it could give an interesting string description. 
For example, due to the additional permutation symmetry Sg, the tachyonic field 
structure can be quite different from the original chain model. The 3-star matrix 
model is also the first one in the class of matrix models having the target space- the 
Dn Dynkin diagramm 
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The gaussian star-matrix model has the partition function: 



\[ dMadMo ex.p[toM^ + UqUq + Mq ^ CaM^ + J2 i^c^^a + UaMa)] (3.1) 



a=l a=l a=l 

The couphng conditions are: 



qPo + 2ioQo + ^to + CaQa = 

a=l 

Pa + '^'taQa + Ua + CaQo = 0, a=l,...q (3.2) 

With the following parametrization of Q-matrices: 

Qo = /+ + aoh + 0.16- (3.3) 

Qa = ba/RaI+ + daIo + Ra^-, a = 1, . . . q 

we arrive at following equations: 

ai = 

2*060; + n + Coi?o = 

2iof^a +Ua + CattQ = 

2to + E ^ = (3-4) 

rto 

2toao + uo + Y^ Cada = 
2ioai + + E Co-^o = 
Solving the coupling conditions we get : 

2? 1 CqWo ri \ 

ai = -^'«o = ^(2^-^ -2«o) 

u Li^^f \ N -2^ 

1 

C?o = ]^(CaWo-2toUo+noE^~<^«X]"|^) 

where A = 4to - E4/*q- 

For quadratic potentials Vo, the basic functions are Hermite polynomials: 



en(Ao) = r?„(Ao) = i/n(^^^7=), (3.5) 

v2cn 

r/o,m(Ao) = i?n(^^7==^),a = 
To calculate the partition function we observe that : 

ho,n = ha,n = Ra (3-6) 
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giving: 



Z = const l^^flj^j ' (3-7) 

const in our case is the exponent exp[— l/4(j4aQ + J2a'^a/ta)] obtained by shifting 
the matrices M„ so that the Uniar terms in the potential vanish. 

The partition function and the first simplest 1-point correlation functions as 
TrQ^jTrQa can be calculated by direct integration of the original integral (|3.1|). 
This is not the case for more complicated 1-point correlation functions as those given 
by (2. 14). Instead with the Q-matrix approach this is easy , using the explicit form 
of Q-matrices.In the Dyson-Schwinger approach ,this calculation is also possible but 
only on the sphere. 

As an example we give the result for the 1-point correlation function (when uq = 
Ua = and gq = da = 0): 

n m n-2km-2l n'm'(-l)'^2"('+*^) 

Tr«Q3- = Z E E E iy,,«,(„_2.-i).{™-2i-i). 

In the large limit we have l + k = Oovl = 0,k = 0,i+j = ""^"^ and the previous 
formula simplifies to: 

For m = and n = 2r we obtain the correlation function in genus 0: 

= , iV+^(6.r (3.10) 

r!(r -|- Ij! 



We have sum over i = r, j = 0. 



4 The gaussian model on Bethe tree 

Kazakov and Migdal Q have obtained the so-called induced "induced QCD"- a ma- 
trix model embedded in the regular D-dimensional lattice. For gaussian potential , 
the model was solved by Gross |^]. In the limit N ^ oo the Kazakov-Migdal model 
with generic potential is equivalent to the matrix model with a Bethe tree target 
space 

Because the model was studied only with the saddle-point method, it is interesting 
to study it in a different framework , that of Q-matrices approach. It gives higher 
accuracy in studying different regions of Bethe tree and also permits computations 
in higher genera. 

Our model is the inhomogenous version of matrix model on Bethe tree, at every 
level of branch we assigne a specific partition function and propagator. 
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The gaussian matrix model on Bethe tree is {i = 2j): 

Z = JdMi exp Tr[^(iiM2 + mMi) + CiMiMj] 



(4.1) 



<ij> 



where < ij > denotes the permitted Hnks of Bethe tree, and Cj are the couphng 
constants of the i-th level branch. 
The coupling conditions are: 



qPi + 2tiQi + Ui + CiQi-i + {q - l)ci+iQi+i = 
qPi+i + 2ti+iQi+i + Ui+i + Ci+iQi + {q- l)ci+2Qi+i = 

Introducing Q-matrices with the form: 

Qi = Sil+ + ttilo + fiC- 
Qi+i = bi/RiI+ + dilo + RiE-, 

we get from the coupling conditions the following equations for coefficients: 

2tiSi + ci^ + {q- = 

2tiai + Ui + Cidi-i + (g - l)ci+idi = 



(4.2) 
(4.3) 



2tifi + + CiRi-i + {q- l)ci+iRi 







(4.4) 



2ti+iRi + Cj+i/j + {q- l)cj+2/i+i = 
2ti+idi + Ui+i + Q+iOj + {q- l)ci+2ai+i = 

2U+i^ + ^ + Ci+iSi + {q- l)ci+2Si+i = 

Let suppose that all coefficients for various i are equal Cj = c,ti = t,Ui = u.ln this 
case the equations reduce to the following set of equations: 



f{Si) + g{Ri-i) = 0, f{Ri) - g{Si) = 
f{ai)+K = 0, f{di)+K = 



(4.5) 



/( 



Ri 



2qnt 1 
~^R'i 



= 0, /(/.)- ^1=0 



where the constant K = (q — and the functions f,g are: 

2' 



Xi-l + 



2{q - 1) - 4 



Xi + {q- 1) Xi+i 



C Xi Xi-\.\ 



(4.6) 



4.1 Fractal regime 

In the scaling limit Ri R, Si 
n = l): 



S,bi 



b, fi f,di ^ d,ai ^ a and (we take 



RS 



4t^-{qcy' 
7 



qc - 2t 

a = a 



4*2 - (gc)2 

This case represents the hmiting case of a matrix model on the fractal curve of 
Bethe tree. We can define the free energy on the unit length: 

F;_ = logi?S = log^^^^^ (4.8) 

We observe a singularity of the free energy at (2t/c)^ = 1. This critical point is 
the analog of critical point for c = 1 matrix model (or g-multimatrix model) at the 
self-dual radius = 1 We must consider the physical domain when 2t < qc {t < 
to have a well-defined path-integral). We choose c < 0. We see that the other region 
2t > qc, c > is not reached. 

We can calculate the 1-point correlation function (for simplicity we take u = 
and genus 0) 

(tVM2*^\ = (TrM/^i) = T¥(Q2^-) = 6'^Tr(/+ + e^f^ 



\ ' /o \ '"^^/o 



(2A;)!iV^+i / -qt 



A;!(A: + 1)! \At'^ - {qc^ 



(4.9) 



4.2 Asymptotic regime 



We have considered the fractal curve or a tiny strip of surface which is filling densily 
the extremity of the Bethe tree. We can study a larger strip which tends asymptoti- 
cally to the fractal curve. 

In the large limit we can scale the coefficients Ri — > ^/N Ri(x),Si \fN Si{x), bi 
Nbi{x),fi — > ^/Nfi{x),di — > di{x),ai — > ai{x), where x = n/N. We can see now that 
all second terms in the equations ( [4. 61 ) are proportional with x, < x < 1 and can 
be considered as perturbations. We neglect the function g. In this case we can solve 
the recursion relations ( [4.6| ) with the result: 

R 'i — ^ J Si = s , 

a, = r'-^K,d, = s'-^K, K=^-u (4.10) 
q^c q'^c It 

hi = b'r\fi = f 

where r = r_|_, s = r_ or viceversa with r± being the solution of second order eq.: 

(g-l)V2 +(2(g-l)-— )r± + l = (4.11) 
and b = f = b- or viceversa with b± being the solution of second order eq.: 

{q - + {2{q - 1) + ^ - ^)6± + 1 = (4.12) 

We see that the alternating Qi,Qi+i in the asymptotic regime can be inter- 
changed; hence it does not matter if i is odd or even. 
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We can define also in this regime the free energy on the unit length: 



F,,,;+2 = log((7 - irRiSi = (4.13) 

(q must be bigger than 2 to have at least one branch). The free energy is and is 
different from Fj^ac- This is understandable because -F/rac is proportional with x, 
but we have considered the case when x = 0, hence -Fj,i+2 = 0. To see if -Fi,i+2 really 
tends to Ff^-ac we must include the perturbation in x. 

As we said before, our model differs from the one studied by Gross and Boulatov 
[^1^]. Their model is homogenous ;for the Bethe lattice with coordonation number 
2D their partition function is: 

Z = J dMiexpTr[-^^M2+ ^ MiMj] = J dMZ{Mf^e^p I-^-TtmA 

i <ij> ^ ^ 

where the partition function of a branch Z{M) satisfies the equation: 

Z{M)= J dM'Z(M')^^~^ exp Tr(- + M'M) 

Our model is inhomogenous. Every branch of different level i has a different partition 
function Zi{M). Hence the total partition function is: 

Z = j dMZ{M)''exp (^-^TrM^^ 

and Zi{M) satisfies the equation: 

Zi{M)= J dM'Zi+i{M'y-^Zi_i{M')expTr{--M'^ + M'M) (4.14) 

We see that the partition function for i-th level branch is expressed not only in terms 
of higher level branches (i + l), but we have also the back-reaction on the lower level 



branches {i — l).We also observe that the eq.( |4.14 ) is different for i odd or even. 



Solving the equation is equivalent with the first two equations ( ^.61 ). Solving 
them as recursion equations we get the power-like solution (4.11) Zj{M) = r^"^^^^ , 



where r satisfies the second order equation ( [4.111 ). If we solve ( |4.6| ) as differential 
equations (after proper scaling when |(i + l) — 2| <^ i) we get the exponential solution 
Zj[M) = exp(jrTrM^) where r is (from equation ( |4.11| )): 



The signs alternate for i odd and even. This result can be compared with that of 
homogenous model if we identify t/c = m'^/2,q = 2D. The partition function per 
branch for homogenous model behaves as Z(M) = exp(— aTVM^), where a is: 

m'^{D - 1) ± DJm* - A(2D - 1) 

ct± = 7^ 

2D -I 

We have an interesting property in the asymptotic regime:the coefficient Rk near the 
point ^^(i/c)^ = 1/2 has a slow oscilation with the period T^fc- 
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The coefficient Rj. is directly related with the free energy of branches of level 
between k and k + 1: 

Fk,k+i = log(g - if'^Rk 

Sk has a complementary oscilation such that Ffc^fc+i(i?fc)+Ffc+i^fc_|_2(S'fc) = ^^^^+2 = 0. 
This behaviour is typical to an antiferromagnet: Rk,Sk are like spin-up and spin- 
down configurations which group pair-wise having a total energy zero. 
We take for convenience i?o = 1, -Ri = 0. The coefficient Rk is: 



Rk 



where r± are the solutions of the equation ( [4.11|) . We introduce the notation /? 
^-j-(t/c)^ and consider the region < /? < 1. In this case: 



r± = J^e^-, arctan. = ^^^ni (4.15) 

In this region R^ = {q — has a fast decaying in amplitude (for q > 2) and 

an oscilatory character with pulsation u> . For /? ~ ^ the expression ( 5.16I ) is singular 
and the pulsation is u ~ ^. If we take /3 = ^ + e then co = ^ + Auj and Auj ~ e. This 
property induces a modulation of the oscilation with the period : 

For (3 1/2 the modulation dissapears because the period T^^ 00. 



5 The q-Potts-like model 

q-state Potts spins are an interesting generalization of the Ising model {q = 2) . On 
planar random lattice they were studied first time by Kazakov [Q] . The cases q ^ 
and g — > 1 represent the models of tree-polymers, respectively that of percolation. 
The q = 2 case , that of Ising model on random lattice, can be expressed in terms of 
the 2-matrix model. 

The partition function of q-state Potts model on a random lattice is: 

n n n 

Zig,l3,H)=J2g- ^ ^ exp[-| J] Gi;)(5.,.^ - 1) + - 1)] (5.1) 

" {G(")}{(t} k,j k 

where the summ run over all triangulations with n triangles {G*-"-*} and all spin 
configurations {cj}. G^"^- is the adjancency matrix of planar lattice with n vertices, 
[3- inverse temperature, H- magnetic field . 

For zero magnetic field, the partition function can be expressed in terms of the 
matrix model: 

^ g <? g 

Z= Y[ dM^ exp[2c J2 MaMp - J2 iK + g/^Ml)] (5.2) 

a=l a>l3 a=l 
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The partition functions (|5.2D and are equal due to the equivalence of the 

Feynman graphs of matrix model on dual lattice with inverse temperature f3duai = 
log(l + g(e^ — 1)~^) and the Boltzmann weights of Potts model on the original lattice. 
Introducing a new matrix Mq the previous integral can be rewritten as: 

Z = y n dM^dMQ exp[cMo ^ - Ml/2 - ^ {Ml/2 - ~g/?,Ml)] (5.3) 

«=1 a=l 0=1 

with (? = c/(l + c),^^ = 5^/(2(1 + c))^. The coupling constant c is connected with 
the inverse temperature (3 by the formula 

We consider the star-matrix model with partition function: 

^9 q 
Z= Yl dMc^dMo exp[toMl + uqMo + Mq ^ c„M„ + 
0=1 0=1 
q 

+ ^ (s^M^ + taMl + n„M«)] (5.4) 

a=l 

The coupling conditions are: 

q 

qPo + 2tQQo + no + ^ CaQa = 
Pa + 'iSaQa + 2taQa+Ua + CaQo = 0, Ct = 1, . . . Q 

With the following parametrization of Q-matrices: 

Qo = 1+ + o-qIq + ai-^- + 02-^^-2 

Qa = ba/RaI+ + + Ral-,a = I, . . . q 

we arrive at following equations: 

^SaRa + Caa2 = 

oi = 

+ n + Cai?Q, = 
•iSaida + 26q,) + 2tada + Ua + CattQ = 




2to + E 



2toao + no + E Cad« = 
2ioai + qn + Y, CaRa = 



(5.6) 



5.1 Symmetric case 

We solve the special symmetric case when Sa = s,ta = t, Ua = u,Ca = c. In this case 
Ra = R, da = d^b/y. = b . 

We can express all the coefficients in terms of only two of them R and d: 
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aa = i?',ai = - 

c 2to 

uo + qcd 2toR 

ao = ,b = (5.7) 

2to qc 



R and d satisfy a system of 2 non-linear equations: 

t cq . n 



f? - 9*^-72 I / get (9c)^c gCT (gc)cno 



Hence all coeficients of the Q-matrices can be expressed in terms of the R^ = R 
coeficient which satisfies a third-order equation: 

R' + R'r^. i ^-^^^^^^ + 6.(cuo - 2tou)] - = (5.9) 

(12sto)2 \ 4to ^ ;y 4to(,l2sto)^ ^ ' 

For symplicity we choose 6s = 2to = 2t = —I^uq = ^,c^ = 1/q. If we denote 
z = R/qc we have instead ( |5.9| ) the equation: 

z3-z2n-M/2 = (5.10) 

Two roots coincide zi = Z2 = -z* = 1 at the critical point when the " cosmological 
constant" = 3/2. Near critical point, the variable R related with the free energy 
will scale as i? — i?* (n — li*)^/^^^'^^^) for the {p,q) matter models coupled with the 
2d gravity. 

Expanding u and z near critical point: 

u = u^, + z = z^, + Zb^ 

we get in the lowest order of b the relation /i = SZ^. This means that the variable R 
scales s& R — R^ ^ (u — n^,)^/2.This critical point (the continuum limit) corresponds 
to the pure gravity model or 1-matrix model. The results remain true for arbitrary 
s, ^) ^*0) c in the symmetric case. 



5.2 Non-symmetric case 

We can write the system of equations ( |5.5[ ) , (^^) in a different way which will permit 
to remain with a single type of variables X^'- 

Xa = 3Sada+t'^ (5.11) 

We denote by: 

2*0'"^"^" 4io'""" 4io 
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Then we can express the variables Ra in terms of Xq. from the system: 

XaRa + Y] CafiRp = a = l...q (5.12) 



and also the variables ha in terms of Xa- 

2baXa + n 2b 13 



p +^c,/3^ = 0, a = l...(z (5.13) 



We remain with the system: 

^ + 3sai2ba) + < - ^ + ^ ^(X^ - t;^) = 0, a = 1 . . . g (5.14) 



We have the suplimentary constraint which can be imposed on the variables X^: 

^^R^ = const, a = 1 . . .q (5.15) 



Co 



When we tend to the symmetric case with Ca/s = c,t'a = 0,6sa = 2to = — 1,Mq 
u + cP{q — l)^ we get from the system ( ^.15 ) the following equation: 

^<-^' + ''<^-"''- 2(X + cfa-l)) -'' = ° 



With the notation z = R/q = —1/{2X + 2c{q — 1)) we obtain the equation (5.10). 
q=2 case 

We argue that the critical behaviour of the case q = 2 coincides with the Ising 
model on random cp^ lattice. 

Solving the system ( ^.13[) we get: 

n{X2 — cr) n{Xi—c/r) 



where: 



X1X2 - X1X2 - 



Rl {ClS2V'^ 



R2 \C2S1 



(5.16) 



Because i?i,i?2 are: 



{X2-C) 2n (Xi-c) 2n 

XiX2-c2 4to' ' X,X2-C^4to 

the relation ( p. 16 ) can be rewritten as: 

X2-C 

r = 

Xi-c 

We remain with the system: 

X? 2c X2-cr , t? t'n 

X2-3sin—j- = ^ + ^.17 

3si 3S2 X1X2 — C^ 6S\ 6S2 

Xl 2c Xi - c/r , t'^ t\ 

6S2 0S2 X1X2 — C^ 6S2 oS\ 
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We point out the great similarity between this system and that of the Ising model 
on (f)^ lattice [ p!o[] (which corresponds to the case q = 2 Potts model). To show this, 
we integrate the intermediate matrix Mq in the relation (^!J).We get the two-matrix 
model: 

r ^ 
Z= dMidM2 exp[^ {saM^ + t'^M^ + u'^M^) + 2CM1M2] 

with the previous notations for t'^,u'^,c = Ci2- 

Solving the coupling conditions we remain with three equations: 

X,X2 = c' + —, 

^ + (^ ^)X2 + u[ = -^ + 2c-^, 5.18 

6si 6S2 2c 6si 6s2 

Xl ,2c 3s24i?,,^ , if t\ 

— + — )Xi +4 = ^ + 2c^ 

3s2 3s2 2c 3s2 3si 

Introducing the expression of 2R = ncj (X\X2 — c^) in the last two equations ,we see 
that they differ from the equations ( ^.18 ) only by the terms containing the r quantity. 



We expect that these terms are only an artefact of the different basis of orthogonal 
polynomials we have chosen and that they do not change the critical behaviour of 
the free energy. The variable R scales at the critical point as i? — i?* ~ (u — u^^l"^ 
(we take n'^ = n, n2 = 0). 
q=3 case 

This is the first non-trivial member of the q-Potts-type set of models. It can not 
be derived from the two-matrix model or other more complicated chain models. It 
represents the matrix model on the Dynkin diagramm D-^. 

We have not managed to derive the critical behaviour. However, for further devel- 
opments, we write down the system of equations which gives the critical scaling. 

We introduce the function: 

^123(1,^,^2) = X2X3 - C23 +ri(ci3C23 - C12X3) +r2(ci2C23 - C13X2) (5.19) 

The indices 123 of the function 1^123 are related with the indices of the variables 
^o,Cq,/3. The function is symmetric only in the last two indices: ^123(1, '''1, ''^2) = 
^132(1, r2,ri). 

We also introduce: 

Y = X1X2X3 + 2C12C13C23 - (0^3X2 + cigXi + cj^Xs) 



Then the coupling conditions ( 5.14| ) in this case are (for simplicity we consider 
t'a = 0): 

^1 J. 2ci2 y ^ 2ci3 y „ ^123(1, ri,r2) , 

IT- + — -^2 + - — X3 - Ssm +ui = 0, 

3.S1 3s2 3s3 Y 

Xl 2C12 2C23 Y o ^23l(l,r2Al,l/n) ^ , n on^ 

IT- + ^ — ^1 + ^5 — ^3 - 3sin + U2 = 0, 5.20 

3s2 osi 3s3 Y 

xl ^ 2ci3 ^ ^ 2C23 ^ „ ^312(1, ri/r2,l/r2) , 

— + - — Xl + - — X2 - Ssm + = 0, 

3s3 3s3 3s2 Y 
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The variables ri,r2 can be written in terms of X^'- 



^23(1. 1,1) t'l23(l. 1, 1) 



6 Conclusions 

We have studied the inhomogenous matrix model on the Bethe tree and have ob- 
tained similar results with the homogenous model in the saddle-point method. We 
have two regimes: the fractal and the asymptotic .In the asymptotic regime, we get 
the exponential behaviour for the partial partition function Zj[M) = exp(jrTrM'^), 
where r satisfies a second order equation like in the homogenous case. For large j- the 
level of the branch- we expect that the properties of the model become independent 
of j and is a transition from the asymptotic regime to the fractal regime. Also in the 
asymptotic regime when /? = = 1/2 we have a slow oscilation of the free 

energy with the period T ~ (1/2 — (3)^^. 

For the q-Potts-like model with arbitrary q we write down the general coupling 
conditions. For the special cases q = 2 and q = 3 we solve the coupling conditions in 
terms of only one type of variables Xa. For q = 2 we have argued that the system 
has the critical behaviour of the Ising model on cp^ lattice . 
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